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One of the main difficulties in studying Quantum Field Theory, in the perturbative regime, is the 
calculation of D-dimensional Feynman integrals. In general, one introduces the so-called Feynman 
parameters and associated with them the cumbersome parametric integrals. Solving these integrals 
beyond the one-loop level can be a difficult task. Negative dimensional integration method (NDIM) 
is a technique whereby such problem is dramatically reduced. In this work we present the calculation 
of two-loop integrals in three diferent cases: scalar ones with three diferent masses, massless with 
arbitrary tensor rank, with N-insertions of a 2-loop diagram. 
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Perturbative calculations in Quantum Field Theory are often a hard task, specially if one does not have a suitable 
approach to tackle the problem. Among the several techniques available in the market the most popular is the 
Feynman parametrization jj]; in fact, if one is clever enough hefty calculations (at 2-loop level) can be done 0. 
However, in our point of view, this is not the most adequate nor elegant method to solve Feynman integrals whether 
■ one is considering covariant or non-covariant gauges. 

On the other hand, Chetyrkin et al developed the integration by parts in configuration space (associated with 
Gegenbauer polynomials) and performed (at 4-loops) even heftier calculations ||. However, their technique has a 
drawback: if the diagrams have more than two external legs manipulation of Gegenbauer polynomials becomes very 
difficult H to handle. 

Mellin-Barnes' contour integration is a third option in the market. Each propagator is Mellin-transformed B 
a very simple step - and using Barnes' lemmas and summing over the residues, it is possible to write the Feynman 
integrals as hypergeometric functions or hypergeometric-like series. Smirnov |H) solved the massless double-box with 
Mellin-Barnes approach. Such Mellin integrals are parametric-like integrals - though they are much simpler to solve 
than Feynman-like ones because Cauchy theorem can be applied straightforwardly. 

NDIM H is a technique where the parametric integrals do not appear. We start up with a Gaussian integral, 
which is well-behaved, perform Taylor expansion and solve systems of algebraic equations. All the calculations can be 
done analytically and the results are given for arbitrary exponents of propagators and space-time dimension D, as in 
the standard dimensional regularization |?]]. Even integrals pertaining to the trickier non-covariant gauges, such the 
light-cone |8 and Coulomb j^] ones, can be performed using the same approach (however, we will deal with them in 
another work). In the following sections we will show how this can be done. 

The outline for our paper is as follows: in section II we consider the simplest two-loop diagram, our workhorse, as a 
pedagogical example and to present, in a clean way, the technique of negative-dimensional integration. Section III is 
devoted to the referred diagram but now with tensorial structure. The NDIM can handle vector, second rank tensor 
and higher integrals, all at the same time. In the fourth section we replace massless propagators by massive ones and 
in section V, we consider the massless diagram with N-insertions of the same type. In the last section, VI, we present 
our concluding remarks. 

II. SIMPLEST TWO-LOOP DIAGRAM. 

To make things clear we begin with the diagram of figure 1. In a massless scalar theory it is represented by, 

J (q 2 W)(p-q-kr- { 1 
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Negative-dimensional approach. 
Consider the gaussian integral, 



G = J d D q d D k exp [~aq 2 - (3k 2 -<y(p-q- k) 2 ] , (2) 

where (a, (3, 7) are such that G is well-behaved. We will see that it is the generating functional of negative-dimensional 
integrals, see eq.(|^). Integrating over q and k is very easy, 

G= (*I) D,2 ap (-2Ei p A (3) 



where A = a(3 + «7 + /?7. Expanding (^) in Taylor series and using a multinomial expansion for A we get, 

_°° n ni23 «ni24_,ni34/' T) 2^ni 

^ ni!n 2 !n 3 !n4! 

ni,n2,n3,n4— 

where due to our multinomial expansion, 0,234 = — n\ — D /2, and we define 

n 12 =n 1 +n 2 , n 123 = n x + n 2 + n 3 , 

and so forth. 

On the other hand, Taylor expanding ^ , 

G = £ ( " 1) ^y^ 7 ' / d D q d D k ( q 2 T(k 2 y( P q kf\ (5) 
i,j,i=o ■ J 

one generates the negative-!? integral. Now, comparing (Q) and (^) we solve for the integral above, 

A(i,j,l) = J d D qd D k (q 2 T(k 2 y(p-q-k) 21 (6) 

* D iW V> (P 2 )" 1 (-"1234 -0/2)1 , 

~ (-1)1+3+1 m'no'm'n^' ° ra «3,i°ni24 ,J°n 13i ,lO n23i ,- ni -D/2, [I) 

n 1 ,n 2 ,n 3 ,n i =0 L ' A ' J ' 4 ' 

where the Kronecker deltas give rise to a system of 4 equations and 4 "unknowns" . Plugging the solution into (Q) 
provides us the result of A(i, j, I), namely, 

Ad 7 1) = ^r(i + pr(i + jjrq + pr(i - a - g/g 

T(^-^-D/2)^(l-j-D/2)T(^-l-D/2)T(^ + ay P,, [) 

where a = i + j + l + D. However, the above result is valid on the negative-dimensional region and positive exponents 
of propagators (i, j, I > 0). To bring it to our real physical world we must invoke the principle of analytic continuation. 
This is a quite simple step Jl0[ : group the gamma functions into Pochhammer symbols and use one of its properties, 

Doing so, we get the result in our positive-dimensional world, 

A AC (i,j, I) = <K D {p 2 Y{-i\i +j + D/2)(-j\j + k + D/2)(-k\i + k + D/2)(p + D/2\ -2a- D/2), (10) 

and negative exponents of propagators < in Euclidean space. Observe that the result is symmetric in the 

propagators exponents reflecting the symmetry of figure 1. 

The recipe for calculating Feynman integrals using NDIM technology is quite simple: i) to each loop write a 
Gaussian integral whose arguments are the propagators of the diagram in question; ii) complete the square(s) and 
integrate; iii) take the original Gaussian integral, Taylor expand the exponential and change the order ^ <-+ J ; this 
operation generates negative-dimensional integrals; iv) the equality of these two expressions must hold, so one can 
solve for the negative-dimensional integral and gets n-fold series involving Kronecker deltas; v) such Kronecker deltas 
give rise to a system of linear algebraic equations, in most cases it has not a unique solution, since it is a rectangular 
matrix p!l| ; vi) Plugging the solution(s) into that series representation provides the result of negative-dimensional 
integral - sometimes in massless cases there appears degenerate solutions |L0| ; vii) analytically continue the referred 
result to positive-dimensional region using the above property of Pochhammer symbols. The whole procedure is quite 
simple, and we will show that in cases of interest. 
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III. TENSORIAL STRUCTURE. 



The previous result was obtained with amazing easiness: arbitrary negative exponents of propagators, positive 
dimension and no numerical calculations. The reader can rightfully ask: Does NDIM work for tensorial numerators 
as well? The answer is yes Jl^] . We need to modify only one thing. 

Consider the integral, 



B(i,j,l,m)= I d D qd D k (q 2 y(k 2 y(p- q -k) 2l (2q-pr. (11) 



From, 



G T = J d D q d D k exp [-aq 2 - (3k 2 - 7 (p - ? - kf - 2<f>q ■ p] (12) 

71-2 \ ___ ( a h + 2 ^ ~ P<P ~ l<t> % 2 



a; cxp l a p 

(_1 \«120"2 ri ,ni57 flni2356^,«124670" 2+2r * 34 f' T1 2\ni234 

< D E — ,,,,,,, — (-^34 -13/2)!, (13) 

m ,ri7— 



i,j 7 l,m—0 



d D q;d D k (q 2 y(k 2 y(p- q ~k) 2l (2q-p) m . (14) 



As we did in the previous section, solving for B, leads to 

N _ 77^/l.yl/lml ~ 2 "'(-l)""(p 2 )""34(_ Tll2 3 4 _£ ) /2)! jf 
l ,J " ] " (_i)<+i+i+m 2^ =q mlnalnaln^lnBlnalnr! ^r.i'We^n^,* W 

x ^n„ 2+ 2„3 4 ,m^ri 56 7,-n 12 34-r>/2- 

Observe that now the system generated by the Kronecker deltas does not have a unique solution, since there are 
seven "unknowns" and five equations and, by this very reason, two of them will be left undetermined. There are 
(Cj = 7!/5!2! = 21) distinct ways of solving this 5x7 system, but 5 of them have no solution at all. 

In a previous work |Q we did show that all non-trivial solutions are legitimate and lead to the same result for the 
Fcynman integral in question. Here the same occurs - but we will not prove it. 

Note that the tensorial sector of the Feynman integral B(i, j,l,m) is contained in the factor (2q ■ p) m and for this 
very reason the exponent m can not be analytically continued to allow for negative values. In other words, we must 
invoke the principle of analytic continuation for three exponents of propagators leaving the fourth, m, untouched. 

All the solutions will give rise to a double series of hypergeometric type, since we have a 7-fold series and only 
5 Kronecker deltas in (|l5|). However, from the theory of hypergeometric series |l3|] we know that when one of its 
numerator parameters is a negative integer, say — m, the series is truncated and has only m terms. For this reason 
we will consider a solution that is obtained when {713, 714} are left undetermined. 

nr - i \ (-m/2|rc 34 )(l/2 - m/2\n 3i )(D/2 + jMQP/2 + l\n 3 ) 

B{i,3,l,m) = g B > 7— — ; ; 77- : =777^ : — : — : , (16) 

^— ' _ (1 + <J - m 7134HI - 1 - m - D A n34jn 3 !n4! 

713,714— x 



where 



(-n) D (p 2 r'2 m T(l + i)T(l + j)r(l + Z)r(l -a'- D/2) 
9B r(l-j-D/2)T(l-l-D/2)T(l-i-m-D/2)T(l + cr' -m) 1 { ' 

and we use the relation Jl3l, 



(a\2b) = 2 2b (a/2|6)(l/2 + a/2\b). 

Note that for positive to, which is the relevant condition, the series (|6|) is always truncated (when it is even the 
first factor in the numerator is the responsible for that, whereas for to odd the second factor truncates the series.) 
Analytically continuation of <?e gives, 
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-AC 



n D (pY2 m (-i\ -3-1- D)(-j\j + I + D/2)(-l\j + I + D/2)(a' + D/2\i + m- a'), 



(18) 



and the final result, in positive dimension, is given by the series in equation (16) times g^ c , 

(-m/2|n 34 )(l/2 - m/2\n M )(D/2 + j\n 4 )(D/2 + l\n 3 ) 



B AC (i,j,l,m) = g§ c J2 



ri3,n4— 



(1 + u' — m|n-34)(l — % — m — D /2\n^)n^-ri^\ 



(19) 



Observe that for m = we obtain the scalar case (|To|) , for m — 1 an integral with vector numerator, for m = 2 second 
rank tensor and so forth. The results, of course, are contracted with the external momentum p+ The astonishing 
point is that all these new results are contained in the same formula, namely equation (fL9h . 



IV. MASSIVE PROPAGATORS. 



NDIM is a powerful technique. It gives, simultaneously, vector, second rank tensor and higher order integrals. A 
second question one could ask is: Does NDIM work for massive propagators as well? The answer is also yes and we 
need to do only slight modifications. 

Let our generating function, corresponding to diagram of figure 1 where now the virtual particles have distinct 
masses, be 



G m = J d D q d D r exp {-a(q 2 - m\) - f3(r 2 - m\) - 7 [(p - q - rf - mf\ } (20) 
i\j\k\ 



where, 



M(i,j,k)= I d D qd D r (q 2 -mlY(r 2 -m 2 2 y[(p-q-r) 2 -ml] k , (21 



using (H) we get, 



G m = exp (am\ + (3m\ + -ymfjG, (22) 
and following the general procedure, as in the previous cases, one can write the integral as, 

kau nDlljlkl (-»4 - D/2)\{m 2 )^ (m|g K)" 3 (-P 2 )" 4 , x x x r „s 

JVl[l,J, K) — \i + j + k Z^i 1 1 ni456,i n2457,i n3467,fc n 4 567,-£'/2- l^Jj 

^ L> ni ,...,n 7 =a ni....n 7 . 

In this case, the Kronecker deltas give rise to a 4 x 7 system of linear algebraic equations. We have 35 possible solutions 
for such system but 15 of them have no solution at all. So, we are left with 20 triple series, following the prescription 
of summing the ones which have the same variables |ll],|8| (it is equivalent to say that we sum the ones which have the 
same region of convergence H), we get four possible triple series (of hypergeometric type) representing the Feynman 
integral M(i,j, k), 




For the first one we have eight solutions, the second, third and fourth have four, so we have 20 possible solutions of 
the system generated by the Kronecker deltas (8 + 4 + 4 + 4 = 20). Since the last three have the same form, due to 
the symmetry of the diagram in question, we will consider only one of them; the others can be obtained changing 
masses and exponents of propagators. 

We quote only the results, where the analytic continuation process has been already carried out. The first triple 
series is, 
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Mi(i,j, k, {z}) = /i^ 3) (-fc, 1 - k - D/2; 1 + i + D/2, l+j + D/2, l-k- D/2) + (j jfc) + (< <-► k) 

/ 2 ^ 3) (-j - k - D/2, l-j-k-D;l+i + D/2, l-j- D/2, l-k- D/2) + (*«-> j) + (fc «-> i) 



where 



+hT [ c\-a, I- a- D/2; l-i- D/2, l-j- D/2, l-k- D/2), (24) 

h = (-n) D (mly +D / 2 (mly+ D / 2 (p 2 ) k (-i\ - D/2)(-j\ D/2), 
h = {-it) D {m\) l+D/2 {-p 2 y +k+D,2 (-i\ - D/2)(-j\ - k - D/2)(j + k + D\-k- D/2)(-k\2k + D/2), 

h = {-K) D {p 2 T{-i\2i + D/2)(-j\2j + D/2)(-k\2k + D/2)(D/2\ -a- D/2). 

Note that one of the solutions has a factor (0\D/2) = r(D/2)/r(0), so that we have in fact seven terms (in positive 
dimension) instead of eight (in negative dimension.) The second triple series of hypergeometric type is represented 
by a sum of four terms, 

M 2 (i, J, k, {w}) = 3 i^ 3) {-k, D/2; 1 + i + D/2, l + j + D/2, D/2) 

+92F ( c\-j -k- D/2, -j; l-j- D/2, D/2, 1 + i + D/2) 

+g^c ] (-*> -i-k- D/2; l-i- D/2, 1+J + D/2, D/2) (25) 
+g i jf\- ( r, -i-j- D/2; l-i- D/2, l-j- D/2, D/2), 
where the hypergeometric series which appears in the above results 

*#W;7, *,*,{*» = E ^f^ff; . (26) 

ni ,rh2 ,n3— / v i ' 

is a Lauricella function Qj which converges if, 

N<1, and V\xi\ + + VW\< 1) (27) 

and we define, 

* = ( i^ +fc K)' +J/2 K) J+D/2 (^) fc (^»l - D/2)(-j| - D/2), 
.92 = ( i~^ +fc K)^/ 2 (mg)W°/ 2 p/2|i)(-i| - D/2)(-fc| - j D/2), 

M = J^i^TJ+k ( m 2 ) J+ ' D/2 ( m iy +k+D/2 (-J I - D/2)(-fc| - * - D/2){p/2\i), 

54 = 7r D (-m2) CT (-i|i + j + D/2)(—j\i + j + D/2)(D/2\ -i-j- D)(-k\ -i-j-D). 
The variables in J\A± and A4 2 are respectively {x} = {z} and {x} = {w}, where 

2 / 2 2 / 2 2 /2 

z\ = m 1 /p , z 2 = m 2 /p , z 3 = m 3 /p , 

2 / 2 2 / 2 2 /2 

W\ = m l /m z , w 2 — m 2 /m 3 , W3 = p /m 3 . 

Observe that our results, which were obtained simultaneously, agree with Berends' jblj et al ones known in the 
literature (which was obtained by analytic continuation.) In their work, M. 2 (i,j, k, {w}) was calculated first and then 
analytically continued to allow other values of momenta and masses, resulting in M.\(i,j,k,{z\). However, if the 
analytic continuation formula was not known the result in the other region would be difficult to obtain. On the other 
hand, NDIM provides all the results simultaneously, allowing us to obtain also new analytic continuation formulas 

IS 
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V. INSERTIONS. 



A third question the reader could pose is: And for higher number of loops does NDIM accomplish also good results? 
Our main objective is to apply NDIM to higher loops, but up to now we are able to "glue" diagrams with two external 
legs only. 

Consider the diagram of figure 2. It has 4- loops and is represented by, 

C(i,j,k,m,n,s)= f dPq dPr dPl dPt (q 2 ) l (r 2 y (p - q - r) 2k (l 2 ) m (t 2 ) n (r - I - t) 2s , (28) 



in a massless theory. Observe that the integrals in I and t are equal, see eq.(||), to A(m, n, s;p 2 —* r 2 ). Rewriting we 
get, 

C[i,j,k,m,n,s)= [ d D q d D r (q 2 Y(r 2 y (p - q - r) 2k A(m,n, s;p 2 -> r 2 ), (29) 



= Tr D (-m\m + n + D/2)(-n\n + s + D/2){-s\m + s + D/2){(T 1 +D/2\ -2a x - D/2)A(i, j + ci, k), 
= ir 2D (p 2 ) ai+<T2 (-i\i + j + cr 1 + D/2)(-j - ax\j + u x + k + D/2)(-k\i + k + D/2) 
x(ct 2 + D/2\ - 2a 2 - D/2)(-m\m + n + D/2)(-n\n + s + D/2)(-s\m + s + D/2) 
x(ai + D/2\-2a x -D/2), (30) 

where a\ = m + n + s + D, <t 2 = i+ j + k + D. 

Let us study the diagram of figure 3, namely the one that has a single propagator replaced iV-times by the graph 
of figure 1. 

i—N 

c N (v 1 ,v 2 ,...,v m ) = J ...J n^^fe'r 3 - 2 ^)^- 1 [(n-i-qi-n) 2 ]" 3 *, (31) 

i—1 

where Tq = p^ is the external momentum and N is the number of insertions. When N = 1 the above integral reduces 
to A(vi, v 2 .,vz). 

Applying the above procedure we can easily solve such scalar integral, 

Cjf{vx,V2, -,V3n) = A AC (y 1 ,v 2 ,v- i )A AC [v^v^ +a 1 ,v 6 ) x ... x A ac (v 3 n-2, "3N-i + a\ + - + (t N -i,u 3 n), (32) 

where gn = ^3N-2 + v 3 n-i + ^3N + D. 
Water melon diagram. Massless case. 

Recently, the water melon diagram was considered in massive case, whether in four dimensions and analytic result 
(see first reference in fL4| ) or integral representations of it suitable for numerical calculations (see third reference in 

B). 

Applying NDIM we can solve, exactly, such water melon diagrams, see figure 4. Here, we will consider the scalar 
massless case. At two-loop level the referred water melon is the graph of figure 1, our setting Sun diagram. At 
three-loop level we begin with, 



G 



= J d D q d D r d D k cxp [-aq 2 - [3r 2 - -fk 2 - 6{p - q - r - k) 2 } , (33) 
= (£l) D, V-^A (34) 



WAI 



C J V c 

where £ = a(3j + a(38 + a^9 + (3j8. After a little bit of algebra we get a 4 x 4 system which gives 

W 3 (i,j,Wm) = J d D qd D rd D k(q 2 y(r 2 y(k 2 ) l {p-q-r-k) 2m , (35) 

TT ZD / 2 i\j\l\m\T{\ - p 3 - D/2)(p 2 )P 3 



(_l)i+j+!+mr(i - i - £>/2)r(l - j - D/2)r(l - I - D/2)T{\ -m- D/2)r(l + p 3 ) ' 



(36) 



in negative dimension and Euclidean space. We define p 3 = i + j + I + m + 3D/2. Generalizing this result to iV-loops 
is quite easy, 



G 



„ „ i=N 

Wjv(K})= / .../ l[d D q i (q?P(p-q 1 -q 2 -...-q N ) 2v ^, (37) 
n ND /W...VN + i<T(l-p N -D/2)(p2y» 



(-i) s -r(i -v x - d/2) . . . r(i - u N+1 - r>/2)r(i + PN ) ' 

where £1/ = v\ + . . . + I'iv+i and pjv = £f + ND/2. Carrying the analytic continuation out we get, 
W^ C (K}) = k ND/2 {p 2 ) Pn (-pi\2vi + D/2){-v 2 \2v 2 + D/2) . . . {-v N \2v N + D/2)(p N + D/2\ - 2p N - D/2), (39) 

the result for negative exponents of propagators and positive dimension. As far as we know, this result was not known 
for arbitrary exponents of propagators. Observe that when one is allowed to "gluing" diagrams, the result can easily 
be generalized to N-loops. 



VI. CONCLUSION. 



NDIM is a suitable technique to tackle the task of calculating multiloop Feynman integrals. Massless, massive, 
scalar, tensorial, even the ones for non-covariant gauges, such as the light-cone gauge |||| are easily performed. In 
all of them exponents of propagators are left arbitrary as well as the dimension D, just like in plain dimensional 
regularization. Usual parametric integrals are replaced by one Gaussian integral over each momentum flowing in the 
loop, and the main task is to solve systems of linear algebraic equations. Another point that we would like to stress 
is that no numerical calculations are required at all. 
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p — q — k 




Figure 1: Simplest two-loop massless Feynman diagram, setting Sun. 



p — q — k 




Figure 3: 2N-loop massless Feynman diagram, setting Sun with N 
insertions of the setting Sun graph. 




Figure 4: Scalar massless water melon Feynman diagram with iV-loops. 



